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Abstract. The Galactic black hole candidate X-ray binary GX 339—4 spends most of its time in the low/hard 
state, making it an ideal candidate for modeling the assumedly low accretion phase. The radio emission correlates 
very tightly with the X-rays over more than two orders of magnitude in X-ray flux density, suggesting that the 
jet plasma also plays a role at the higher frequencies. We compare the predictions of our jet model, with and 
without acceleration, to thirteen broadband simultaneous or quasi-simultaneous spectra over this changing flux 
history. In addition, we consider a simple standard thin disk which transitions to an optically thin accretion 
flow, in order to account for the assumedly thermal optical data seen in some observations. A solution without 
acceleration cannot describe the data without unrealistic energy requirements, nor explain the non-thermal radio 
spectrum seen during recent radio outbursts. But because of the low disk luminosity, and possibly the assumed 
disk geometry, acceleration in the jet is limited only by synchrotron cooling and can extend easily into the X- 
rays. We present a model which can account for all the broadband spectra included here, by changing only two 
parameters in the jet model: the input power and the location of the first acceleration zone. However, the model 
is most sensitive to changes in the jet power, the varying of which can also account for the slope of the observed 
radio/X-ray correlation analytically. At the highest low/hard state luminosities, the synchrotron self-Compton 
emission from the jet could be detectable with missions such as GLAST, providing a way to test the extent of the 
synchrotron contribution. We conclude that jet synchrotron is a possible way to explain the broadband features 
and this correlation, and discuss ways of incorporating this component into the "standard" corona picture. 
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1. Intioduction 



et al.||l999| — hereafter W99; [Nowak et al] |2002| — hereafter 



The Galactic black hole candidate (BHC) X-ray binary 
(XRB) GX 339—4 displays a wide variety of spectral 
states, and has recently resurfaced at a high flux level 
after being in an extended "off" state since 1999 (see, 



e.g., Kong et al. 2000). It has been extensively, and of- 
ten simultaneously, observed across broad energy bands: 
from the radio (|Hannikainen et al. 1998; Fender et al. 



199£: Corbel ct al. 200C —hereafter COO; Corbel et al. 



optical (Soria et al. 


1999; Shahbaz et al. 


200 1|; Corbel 


fe Fender 


2002— hereafter CF02), to the X-rays (Wilms 
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NWD). These data sets are relatively untainted by con- 
tamination from the companion star. This, and the fact 
that GX 339—4 often transitions between the canonical 
XRB states, makes it an ideal source for comparing the 
emission processes involved in these various states. 



GX 339—4 spends the great majority of its time in the 
low/hard state (LHS) (for a review of BHC states, see 
Nowak 199E), at flux levels varying over roughly three or- 
ders of magnitude in the X-ray waveband (NWD; C02). 
Its LHS spectra are fairly "typical" in that they show a 
flat-to-invcrted radio component, a weak thermal contri- 
bution, and a hard power-law in the soft to hard X-rays 
( ^dziarski etaL| |1998| ; |Fcndcr|pOOl| ; |Rcvnivtsev et al.| pOOl 



NWD). Although no jets have yet been resolved, the radio 
emission shows 2% linear polarization, and the brightness 
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temperature (defined as Si, — QB,j{Tq) where V, is the 
solid angle of the emitting region) requires a spatial ex- 
tension greater than the binary separation (W99; COO). 
At the time of writing, GX 339-4 briefly entered the very 
high state (VHS) and made a transient-like radio out- 
burst which was five times brighter than the brightest 
radio emission ever observed in the LHS. The optically 
thin spectrum had a spectral index a ^ —0.6 {F^, cx j^") 
(R. Fender, S. Corbel, private comm.), which is character- 
istic of a transient bright jet. This makes GX 339—4 the 
first source to show its radio properties clearly in the LHS, 
the high/soft state (HSS) and transient/ VHS. Combined 
with the fact that jets have already been imaged in the 
other three persistent Galactic BH C which share s i milar 
LHS spectra (e.g., 1E1740. 7-2942, [Mirabel et all |l992|; 



GRS 1758-258, 


Rodriguez et al. 


Stirling et al. 


2001 


) , the case for < 



very strong. 

A special feature of GX 339—4 is that its radio emis- 
sion has been observed to correlate tightly with the X-ray 
emission over the entire range of LHS luminosities, and 
also down to a very weak level of emission in the "off" 



state (Hannikaincn et al. 1998; COO; C02). Similarly, dur 



ing changeover to the HSS, the seemingly linked radio and 
hard (20-100 keV) X-ray emission show an anti-correlation 
with the assumed thermal disk emission in the softer (2-12 
keV) X-ray band ( [Fender et al||1999| ; COO). This suggests 



not only that the plasma leading to the synchrotron radi- 
ation plays some role at higher energy, but also that the 



jet is either quenched (Fender et al. 1999) or otherwise un- 



able to form when the thermal disk component dominates 
( [Meier||200l| ; [Merloni fc Fabiaij ^002[ ) . 

On the other hand, the LHS X-ray spectra from 
GX 339—4 show additional features which have been in- 
terpreted as signs of reflection off colder disk material, 
and can be fit by one of several disk-corona models as 
discussed in detail in NWD. In this paper, the authors 
find correlations between the reflection c overing fraction 
and bot h the soft X-ray flux (see also [Zdziarski et al 



1999|; pRcvnivtscv et al.| 2001), and the time lags be- 
tween the soft and hard X-ray variability. Similar to what 
Pottschmidt et al. (2000) report for Cyg X-1 — which a lso 
shows radio/X-ray correlations ( Brocksopp et al. 1999 ) — 
NWD find that these time lags anti-correlate with the 
coronal compactness, and peak shortly after the transi- 
tion from the HSS to the LHS, when the assumed jet 
emission begins to rise. They suggest an interpretation 
that the time lags may be associated with the extended 
synchrotron-emitting plasma. 

While the model fits presented in NWD can well ex- 
plain the X-ray features, and suggest a link between the 
corona and jet, they do not explore the nature of this link 
or how this emission is related to that of the radio. One 
model which can explain the tight radio/X-ray correla- 
tions proposes that the jets directly contribute to the LHS 



X-ray spectrum via synchrotron emission ( Markoff et al 
[200 1|— hereafter MFF ), as occurs in several AGN (e.g., 3C 



Wilson & Yang 2002). However, this model has only been 
applied until now to one source, XTE Jl 118-1-480, which 
may be unusual in that it has extremely weak or absent 
reflection features ( Miller et al.| |200l| ) . 

In this paper, we discuss the application of a jet model 
to the collected simultaneous or quasi-simultaneous data 
sets for GX 339—4 in order to investigate whether it suc- 
ceeds in explaining the observed radio/X-ray correlations. 
We discuss how the model fares against the broadband 
data and compare the model parameters for the various 
data sets to each other, as well as to XTE J1118-I-480. 

2. Data 

The 1981 data set was compiled by CF02, where it is ex- 
plained in fu rther detail. T he X-rays come from the Ariel 6 
instr ument (Ricketts[ 1983 ) , with simultaneou s optical/IR 
from [Pedersenl ( |l98l| ) and [Motch et all ( |l98l| ). 

All radio observations were performed with the 
Australia Telescope Compact Array (ATCA) and, with 
the exception of the 2000 September 15 observation (to 
appear in C02), have been presented previously by COO. 
All observations from COO in which the radio flux density 
was weak (i.e., ^ 1 mJy) were reanalyzed for this paper. 

The X-ray data presented here have been discussed 
previously by NWD and C02. All data were obtained by 
the Rossi X-ray Timing Explorer (RXTE), with the excep- 
tion of the 2000 September 15 data, which were obtained 
by BeppoSAX. 

In this work we are attempting to describe the broad- 
band features of the jet model, and to simultaneously ex- 
plain radio, IR, optical, and X-ray data (taken with a wide 
variety of instruments). Furthermore, we are not yet at- 
tempting to explain detailed features within the X-ray, 
e.g., Fe lines or reflection features. Therefore, rather than 
forward fold the jet models through the detector response 
matrices, we instead compare the models to "unfolded" 
count flux rates. This is not strictly fitting, as we cannot 
obtain chi-squared estimates from unfolded data in this 
way, but for simplicity we will use the word "fit" to de- 
scribe the comparison to the data throughout the paper. 

To unfold the X-ray data, the RXTE data extraction, 
data binning, response matrix generation, and background 
subtraction were carried out in an identical manner as pre- 
sented by NWD. We then fit models comprised of: a multi- 



temperature disk blackbody (e.g., Mitsuda et al. 1984 



66B, [Hardcastle et aljp001| and M87, [MarshaU et al.|pOO^ ; 



Makishima et al, 1986 ) with peak temperature fixed at 
0.25 keV, an exponentially cutoff broken power-law with 
break energy at w lOkeV, and a (potentially) broad gaus- 
sian line with energy fixed to 6.4 keV. (This model is es- 
sentially identical to that discussed by W99). These model 
fits were then unfolded through the detector response ma- 
trix and multiplied by the ratio of the count rate data 
to the fitted model folded through the detector response. 
For the RXTE observations, data was considered only 
for the 3-200 keV range or less for the faintest data sets. 
Furthermore, HEXTE data was normalized to the PCA 
data (e.g., see W99). 
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verse Compton (IC) upscattering by the hot jet electrons 
of both the thermal disk photons, as well as the syn- 
chrotron photons, is also included. 

Our disk model is relatively simplistic, since we are 
currently exploring its limits in terms of contribution and 
there are already many detailed disk models in existence 
(see Poutanen 199^ ; NWD for reviews). We here consider 
one commonly invoked model in which a standard thin. 



optically thick disk (SD) (Shakura & Sunyaev 1973) exists 
down to some transition radius rtr , at which point the flow 
becomes hotter, optically thin and relatively non-radiative 
(e.g., Esin ct al. 1997). This type of non-radiative inner 



"Data from W99, COO, C02 and NWD. All but 00 (which is 
BeppoSAX) are labeled according to RXTE conventions. 
''Date for X-ray observation. Radio taken on 1999.09.01. 
■^Optical data shown in Fig. were taken 3 months prior. 

As discussed by C02, simultaneous RXTE/ASCA ob- 
servations indicate a probable faint background source 
that contaminates the RXTE spectrum — but not the nar- 
row field of view ASCA spectrum — at extremely low flux 
levels. The RXTE observation of 1999 July 29 is assumed 
to be heavily dominated by this contaminating source, and 
this spectrum, multiplied by 0.78, is subtracted, before 
model fitting and data unfolding, from all RXTE obser- 
vations occurring later than the observation of 1999 May 
14. This normalization was determined by comparison of 
spectra obtained simultaneously by RXTE and ASCA on 
two separate occasions (the required normalization con- 
stants for these two observations were 0.73 and 0.83; see 
C02). Note that the 2000 September 15 BeppoSAX obser- 
vation ^s most likely not contaminated by this background 
source, which could easily be half a degree away and still 
affect the RXTE observations, as BeppoSAX has a nar- 
rower field of view than RXTE. A more detailed discussion 
of this background source appears in C02. 

The dates of the observations are listed in Table 1. 
Data sets labeled 20181 were first discussed by W99. Data 
sets 4 0108-02 were first discussed by NWD. Data sets 



disk has als o been invoked fo r low-luminosity AGN such as 
NGC 4258 (lYuan et al.||2002| ). The thermal bump— seen in 
several of the GX 339—4 data sets in the optical band — is 
assumed to originate in the disk, which also provides seed 
photons for jet IC as mentioned above. We assume for now 
that the jet is oriented 90° from the plane of the disk, but 
in the future we will consider other orientations (see, e.g., 



Maccarone||2002| ). 

In their study of the BHC XTE Jl 118-1-480, MFF 
interpreted the extreme-ultraviolet (EUV) upper limits 
(Hynes ct al. 2000) as the hottest, and thus innermost, 
emission from this outer disk. Therefore, they could place 
limits on its maximum temperature and luminosity. They 
derived a location for the transition radius of rtr ~ 10^ — 
10^ 



which is consistent with other models (e.g., Liu et al 



1999| ; |Esin et al.||200l| ). However, no EUV data exist for 
GX 339—4 , and so we constrain rtr with the low-frequency 
X-ray data. We have also made the outer disk slightly 
more realistic than the simple one-temperature compo- 
nent we used before by considering a multi-temperature 
blackbody spectrum as its representation (e.g., Bhakura 
k Sunyaev| |1973| ; [Mitsuda et al.| |1984| ; [Makishinia et al 



1986| ). We then explore its possible geometry using the 
data sets where its signature seems the most prevalent. 

A significant problem for modeling GX 339—4 is that 
most of its system parameters are not well constrained. 
The only report of an orbital peri o d bas ed on optical pho- 
tometry is from Callanan et al. (1992), who found 14.8 



40108- 03, first appearing in C02, were part of the same 
series of observations originally discussed by NWD; how- 
ever, these latter observations were too faint for simultane- 
ous spectral and X-ray variability analyses. Observations 
40105-02-02, 40104-01-08 and 21136001 were part of a se- 
ries of observations to study GX 339—4 in quiescence, and 
are discussed by C02. 

3. Model 

The details of the jet model can be found in MFF, and 
references therein, and we give a brief summary here. We 
consider that an accretion disk may contribute to the near 
infrared (NIR), optical and possibly soft X-ray emission, 
but that the radio through at least far infrared (FIR) are 
due to synchrotron emission from the jets, which could 
possibly extend through the NIR to the X-rays. The in- 



hours . This value remains unconfirmed (though Cowley 
et al. 2002| reconsidered the problem and found a simi 



larly tentative period of 16.8 hours). However, it has been 
adopted over the years to limit the mass function, which 
is also problematic since the companion star is still un- 
classified. Finally, because the system has not yet been 
resolved, the inclination angle is also unknown but limited 
to ^ 60° because of the lack of eclipsing. Based on reason- 
able assumptions, a distance of 4 kpc has been adopted 
(see ^dziarski et al. 1998; W99; but also ^hahbaz et al 



pool 



who argue for > 5.6 kpc) and a mass for the black 
hole ranging from ~ 3 — QMq, where we here use 5Mq. 

We fix the inner temperature of the SD, Tin, based 
on the 1981 data set. We then estimate the outer radius, 
''out, of the SD based on the physical parameters above, 
the best fit inclination angle, 9i, and assuming a com- 
panion of ~ O.4M0, consistent with the 16.8 hour pe- 
riod. The exact luminosity contributed by the SD, L^, is 
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not uniquely determined by the data, and so we can only 
roughly estimate the accretion rate, Md- Because of the 
simplicity of our disk model, this is not a main focus of 
this work, but we want to be sure we are consistent in 
terms of the energy budget. We express in units of 

LEdd = 1-25 X lO^s — 



then similar to, but slightly stronger than, the canonical 



erg s 



For accreting black holes it has been argued that the 
jet power is of order Qj ~ q^M^c^ with an efficiency in- 
ferred to be of order q-^ = lO^'^ — 10^^ (Falcke fc Bicrmann 



199£ ; Meier| |200l| ). We consider Qj a free parameter in our 
model, also expressed in units of Lsdd erg s""'^, and which 
we check for consistency with our limits on M^. While 
the jet formation itself is still a looming question in the 
field, the physics of calculating most of the jet emission 
is relatively straightforward because the flat-to-inverted 
spectrum stems from the part of the jet where it is basi- 
cally undergoing free expansion. We briefly summarize the 
jet component as follows (and see also Falcke & Markofi 
2000| ; MFF). 

At the inner edge of the optically thin accretion flow, 
hot plasma is ejected out from symmetric nozzles, where 
it becomes supersonic. We assume free jets which ac- 
celerate along their axes only due to their longitudi- 
nal pressure gradients. This is the simplest scenario and 
provides a lower limit to the final velocity distribution. 
The velocity field along the jet is thus u niquely deter- 
mined from the Euler equation (see, e.g., Falcke 1996| ), 
and at large distances z from the base has a depen- 
dence 7j(z)/3j(z) oc yj\a{z). Thus, the maximum value 
of the bulk Lorentz factor 7j is roughly determined by 
the distance at which the lowest frequency radio emis- 
sion is emitted and is dependent on other parameters in 
this model but is itself not a fitted parameter. For both 
GX 339-4 and XTE J1118-F480, we find 7j,max 3, which 
is consistent with the low beaming factors suggested by 
the recent work of Gallo et al. (2002). 



We further assume the jets expand sideways with their 
initial proper sound speed, ^sPsC — 0.4c for a hot elec- 
tron/cold proton plasma. Beyond the nozzle region, the 
radius r(z) is related to the distance z by the Mach num- 
ber along the jet, 'jj(3j/%Ps- 

We choose the simplifying assumption that the inter- 
nal energy, dominated by the magnetic field, is equal to 
the bulk kinetic energy of the particles, consistent with 
a magnetic launching mechanism. The plasma is assumed 
to originate in the hot accretion flow and therefore con- 
tains equal numbers of protons and electrons, with the 
exact temperature of the electrons at the base of the jet, 
Te, remaining a fitted parameter. A process such as pair- 
loading via interactions with the disk photon field will not 
be efficient enough in the LHS to be significant (see, e.g., 
Sikora & Madejski ^OOC ). Under these assumptions, the 



magnetic energy density and particle density in each jet 
can thus be defined as B'^{z)/8n ~ 0.25Qj/(c7j/3j7rr(z)^) 
and n^(2) ~ 0.25(5j/(c7j/3jTOpC^7rr(z)^). Taking into ac- 
count the non-constant velocity field, the dependencies 
of the magnetic field, B, and density, n, on distance are 



and 



dependencies for conical jets, respectively 



Johnstoij |1988 



(Plandford fc Konigl| |1979| ; [Hjcllming 
p^'alcke fc Biermami 1995). In this way, the basic physical 
properties governing the emission at each point in the jet 
are fixed after specifying the jet power Qj and nozzle scale 
ro, which determine Bq and ng. 

The thermal particles travel along the jet until they 
encounter an acceleration region which begins at Zgh- An 
acceleration process is inferred because of the observed op- 
tically thin synchrotron power-law seen in the recent tran- 
sient radio events, indicating the presence of non-thermal 
particles. This is a standard process invoked to explain 
many features of AGN jet emission (Marscher & Geai 
1985). We consider the process to be diffusive shock ac- 



celeration, but it could really be one of a few possible sce- 
narios (e.g., stochastic acceleration) which lead to a non- 
thermal power-law with spectral index p (^ oc E~p). 
For standard shock acceleration theory, p ^ 2.0 — 2.5 



(e.g., Jones & Ellison 1991). In order to maintain the non- 
thermal particles out along the jet to account for the entire 
radio spectrum, we assume distributed acceleration after 
Zsh- This is because the cooling timescales after accelera- 
tion are very fast compared to the dynamical timescales, 
a problem which has also required the assumption of dis- 
tributed acceleration in the two best studied AGN jets, 
M87 and 3C 273 ([Meiscnheimer et al.| 1996; Jester et al 
^001| 



but see also Perlman et al, 2001 



We assume that the accelerated distribution extends 
from the peak of the thermal distribution, with Lorentz 
factor 7o ~ 4 • Tg^o- As discussed in MFF, the maximum 
7c occurs when the acceleration rate (assuming Bohm dif- 
fusion) is matched by the cooling rate giving 



7o 



-0.5 ^^sh^ 



(1) 



where Wgh is the shock speed in the plasma frame and B is 
the magnetic field in Gauss. The parameter ^ is the ratio 
between the diffusive scattering mean free path and the 
gyroradius of the particle. It measures how many times 
the particle gyrates per shock crossing, where the most 
efficient acceleration possible is when it crosses the shock 
once per gyroradius, thus giving a lower limit of ^ = 1. 
For quasi-perpendicular shocks ^ < cjePe/ush, although 
in the case of a jet the field is likely quite tangled in which 
case an upper bound on ^ is not known. While we do not 
know the exact value of Ugh, it must be larger than the 
speed of sound at the shock, Cg, and no higher than the 
local bulk velocity. For very relativistic particles ^ can be 
quite large, but conservative values lie in the ~ 10 — 100 
range (e.g., Jokipu 1987). If we define as a reference value 
f = 1^2 100, the maximum synchrotron frequency is 



= 0.29;^^ ^ 1.2 • 10^ 



Hz 



(2) 



where — ^7c max(s-^)/("^cc) is the critical synchrotron 
frequency. This maximum corresponds approximately to 
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the rollover of the power-law cutoff and in order to ex- 
plain the "canonical" 100 keV cutoff for ^ « 100, we 
need Ugh <; Cg, consistent with our physical expectations. 
This cutoff is not dependent on the magnetic field, the 
jet power, or the shock location as long as we are in the 
synchrotron cooling dominated regime. Because we would 
expect XRBs to have similar shock structures, we should 
get roughly similar cutoffs for different sources and accre- 
tion rates. In this sense ^ and Msh provide a reasonable 
physical explanation for the location of the cutoff in much 
the same way the coronal energy fraction and tempera- 



ture determine it in coronal IC models (e.g., Sunyaev & 
Triimp§^^979|). 



The location of the initial acceleration region is de- 
termined by the frequency where the flat, highly self- 
absorbed synchrotron spectrum turns over into the op- 
tically thin power-law produced at the shock. In MFF, 
we showed that from back-extrapolating the X-ray power- 
law, this maximum self-absorption frequency should occur 
somewhere in the IR/optical regime at ~ 10^^ — 10^^ Hz 
for XTE Jl 118-1-480. It is interesting to note that this now 
holds true for several other BHC sources in the LHS in- 
cluding GX 33 9-4 and V404 Cy g (CF02; Brocksopp et al. 
2002, in prep., pallo et alj ^002| , Gallo et al., MNRAS, in 
prep.). 

There are thus seven free parameters relating to the jet 
component, Qj, Tg, Zgh, tq, p and 9i. These are not, how- 
ever, entirely independent of each other. Furthermore, in 
our philosophy of comparing broadband spectral features 
to the jet model, we can choose reasonable constraints for 
a number of these parameters. For instance, if the X-rays 
are due to synchrotron emission, the spectral index of the 



accelerated particles is uniquely determined (ax 



(p-i) 



Starting with reasonable assumptions of the underlying 
disk/jet geometry in order to constrain T^. and tq, Qj and 
9i can then be determined from the radio spectrum be- 
cause they fix its normalization and spectral index, re- 
spectively. The beginning of the acceleration region, Zgh, 
is then uniquely determined by where the turnover must 
occur to link the X-ray to the radio spectrum. For the 
purposes of the spectral models that we describe here, 
the critical fitted parameters are Qj and 6i, as they most 
directly determine the radio spectrum. For 9i, the only 
observational constraint is that it is ^ 60°. 

As described above, there are also additional parame- 
ters describing the cutoff in the X-rays as well as the disk 
component, and we will also argue below for an additional 
irradiation component. However, these are not an integral 
part of our model, in that they do not constrain the condi- 
tions of the radiating plasma in the jet. We include them 
for completeness, and to open the way for more complex, 
combined disk/jet models. 

Some fraction of the emitted X-rays from the jet base 
will either directly interact with, or be scattered by hot 
electrons into, the cold disk, resulting in reflection. This 
calculation is not attempted here, however, since again 
this is very dependent on the disk geometry and we are 



focusing on the role of the jet. This, as well as other issues, 
will be the focus of an upcoming paper. 



4. Results 

4.1. 1981 data set: constraining the GX 339—4 model 
parameters 

Of all the data sets for GX 339—4 considered here, the 
most interesting is that of May 1981 (see Fig. Q), which 
is at the highest flux so far for which simultaneous broad- 
band data have been published in the LHS. The IR/optical 
wavebands seem to indicate the clear presence of both jet 
and disk components, which has not been seen explicitly 
in any other XRB. The IR band shows a first component 
with a negative slope suggesting that we are seeing the ex- 
pected turnover from the optically thick to optically thin 
regimes. The shape of this spectrum at higher frequen- 
cies is hidden under a component which is likely due to 
thermal emission from the SD, as indicated by the sharp 
rise in the optical points. However, if the simultaneously 
measured X-rays are traced back to the IR, they line up 
with the turnover remarkably well, supporting their inter- 
pretation as synchrotron emission. 

While no simultaneous radio observation exists, an es- 
timate of the 8.6 GHz radio flux can be found by extrap- 
olating the radio/X-ray correlation curve to higher fluxes 
(COO; C02 and see Fig. H), and is shown as a point in 
Figs. & §). As mentioned above, several other LHS 
sources show this interesting "turnover coincidence" . In 
the compact cores of most flat-spectrum AGN, th e same 
kind of turnover is observed in the mm range (e.g.. Bloom 
|et al. 1994), but the more compact scales in XRBs push 
this to higher frequencies (e.g., Falcke et al. 2001). 

The clear presence of the optically thick-to-thin 
turnover makes the 1981 data set a strong test of any the- 
oretical model. Furthermore the radio/X-ray correlation, 
illustrated by the full set of simultaneous data shown here, 
suggest that any model for the X-rays from this source 
must address the lower-frequency data as well. This means 
that if — as is usually assumed — the X-rays are due en- 
tirely to IC upscattering off hot electrons in the base of 
the jet/coronal plasma, one has to show that these same 
electrons can reproduce the flat-to-invcrted radio spec- 
trum via their synchrotron emission, or that some other 
very strong coupling exists between "corona" and "jet". 
As Blandford fc Konigl (197£) explain, to fit the flat-to- 
inverted radio spectrum particularly in the cm-band, one 
needs a conical geometry with outflowing plasma, other- 
wise the emitted synchrotron spectrum will be too in- 
verted. We therefore think that a corona which has no 
relation to the jet cannot explain the data for GX 339—4 . 
However, the corona may provide the launching point for 



the jet (Merloni & Fabian 2002) or even comprise the 
base of the jet itself (Fender et al. 1999). Interestingly, 
Beloborodov (199£) found that magnetized plasma mov- 
ing away from the disk with a velocity of 0.3c can ex- 
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plain the weak reflection features in Cyg X-1, which would 
be consistent with this picture. 

In order to understand the IR/optical data, we must 
first consider how the disk may be contributing in this 
region. Unfortunately, addressing the thermal disk com- 
ponent is not a uniquely defined problem for GX 339—4 . 
If we assume the SD component extends up to the first 
X-ray data points and then drops off, we find Tin = 10® 
K, which is the same value which W99 used for their fits 
to this source. We keep Tin fixed to this value for all mod- 
els discussed in this paper. However, a naive fit to the 
rising optical data with a multi-temperature blackbody 
would then require almost the entire TEdd being chan- 
neled into the SD emission. For the LHS, which is defined 
by its assumedly low accretion rate, this seems unrealis- 
tically high. One way around this is to assume a lower 
SD luminosity, and that the optical points are due to the 
irradiation of the outer disk material by intercepted X- 
rays. For the mass and inclination angle we use here, the 
outer disk radius would be on the order of rout ~ 10^ cm. 
We use this as the nominal scale for the X-ray irradiated 
shell, noting that there could easily be a factor of a few 
leeway due to the large uncertainties in the system's bi- 
nary parameters. To fit the optical data, the temperature 
of this single blackbody needs to be Tbb ^ 5 x 10"* K. 
This high temperature is not in conflict with possible disk 
instability models (e.g., Cannizzo et al. 198?) because we 
are observing already in outburst, after the effects of the 
instability would be felt. 

For the models presented in Fig. ^ we fix the SD lu- 
minosity to a characteristic value of Td ~ O.lLEdd, and 
then make up any deficit in the optical region via the irra- 
diation component. The ensuing fraction of the jet power 
contained in the irradiated component then can be deter- 
mined due to the fact that an IR/optical break is observed. 
For the other data sets included in this paper, however, 
this break is not observationally resolved. For those mod- 
els we therefore assume that the luminosity of the irradi- 
ated emission can be linked to the power in the jet, and 
we set Tbb = O.lQj as a canonical value. 

It turns out that the 1981 data set is in fact so luminous 
that it creates a problem for our assumed black hole mass 
of 5M0. Even assuming Td ~ O.lTEdd in fact requires 
m ^ 0.5 (in Eddington units) for the fixed Tin — 10® 
K, which is much larger than one would expect for the 
LHS according to some models (e.g., Esin et al. 1997). We 
could solve this problem by choosing a lower Td and Tin, 
thus having the SD contribution truncate before the X-ray 
data. However, then this only increases the need for more 
intercepted jet power in the irradiated component, which 
we already require to be at the upper bound of physical 
limits. We discuss the parameters in detail for the models 
below, but feel this problem is too unconstrained at this 
time to make any conclusions about the exact nature of 
the disk contribution. 



4.2. 1981 data set, no acceleration 

Assuming that jets are responsible for the radio emission, 
the first issue to explore is how the solution would look 
if the jet plasma never encounters an acceleration region. 
This would be the simplest case model, where a portion 
of the thermal disk or coronal plasma is simply advected 
into the outflow, retaining its thermal or quasi-thermal 
distribution. In this case, the synchrotron spectrum from 
each segment of the jet will have a "hump" shape peaked 
around a particular self-absorption frequency determined 
from the physical parameters for the segment. The com- 
bination of all these humps, each located at a different 
frequency, gives a flat-to- inverted radio to at least IR spec- 
trum which can be seen in Fig. |l]. The slope of the spec- 
trum depends on the amount of expansion and the velocity 
profile along the jet, as well as on 9i. 

As explained in the previous section, the electron tem- 
perature at the base of the jet is a free parameter, and so 
we explore two extremes for this scenario. The first, shown 
in Fig. pL, is at the lowest end of the relativistic scale with 
electrons entering the jet base with initial Tq ^ 2x 10^ K, 
and in Fig. |l|b we show a more relativistic solution, with 
Tc = 2 X 10^" K, the maximum temperature we could ac- 
commodate in our modeling of XTE Jl 118-1-480 (MFF). 
In order to bring the highest frequency "hump" from the 
part of the jet closest to the accretion disk to low enough 
frequencies to come near the turnover in the IR data, we 
must vary the nozzle radius tq and the input power Qj, 
thus controlling the magnetic field Bq and density uq at 
the base of the jet while the temperature remains fixed. 
Similarly we must adjust 9i to fit both the radio and the 
IR together with synchrotron, which changes the power 
requirements because of the beaming. 

With the requirement that the non-accelerated mod- 
els adequately describe the IR data as well as the maxi- 
mum self- absorption turnover, the jet parameters for the 



5 X 10^ 



Qj = 3.7T 



Edd 



fit shown in Fig. ^a are then rg 
and 9i = 8°, the former two being obviously unlikely val- 
ues. The extremely high value of the the nozzle radius is 
necessary to decrease the self- absorption frequency to a 
value low enough to fit the IR data, which because of the 
resulting larger self-absorption also leads to a very peaked 
contribution from the denser nozzle. This does not seem 
feasible under the assumption that the jet originates at the 
inner edges of the accretion fiow, or in the corona above 
it, and the match to the IR data itself is not particularly 
good. The "nozzle bump" is due to our assumption of this 
element of the geometry, as required for jet models of Sgr 
A* a nd LLAGN (e.g., [Falcke fc Markofl^ ^000| ; |Yuan et al. 
^002| ). Without a nozzle region, one could find a solution 



with a somewhat larger 6i, but then this would require 
an even higher Qj because the emission would be beamed 
away from the observer. The value of Qj here is already 
several times higher than TEdd for a 5 Mq central ob- 
ject. The irradiation component can provide the required 
rise in the optical data if the outer disk intercepts a small 
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Fig. 1. Two fits to the GX 339—4 1981 data sets for the case of no particle acceleration in the jet, data references 
contained in CF02. The solid thick line is the total spectrum, the dotted line is the multi-temperature blackbody SD 
disk plus irradiation blackbody contribution , the short-dashed line is the synchrotron emission and the long-dashed 
line is the inverse Compton (IC) upscattered disk and jet photons. In both panels, the IC components do not include 
the outer part of the jet as its contribution is orders of magnitude under the data, (a) The fit for the nearly non- 
relativistic initial electron temperature of 2 x 10^ K and rg = 5 x 10^ r^, (b) The fit for a temperature similar to the 
maximum derived for another XRB source, XTE J1118-F480 (see MFF), of 2 x 10^° K, with tq = 2 x 10^ Vg. 



fraction of the jet power in X-rays, Lbb ~ Q-OiQy Overall, 
however, this model is not so satisfactory. 

The observations are slightly better described, at least 
at low-frequency, for the more relativistic case shown in 
Fig. |l|b. The higher temperature electrons lead to less self 
absorption, giving a broader nozzle bump and requiring 
a somewhat smaller nozzle radius of rg = 1.2 x 10"^ 
with a more reasonable Qj = 1.4 x lO^^iEdd for Oi — 15°. 
This nozzle width is still rather untenable, however, in the 
context of jet/corona models. For this scenario we also 
run into more trouble with the disk modeling, because we 
would require the outer disk to intercept roughly GQj in 
X-ray luminosity to account for the optical data. 

For both cases shown in Fig. |l|, the IC emission in 
the X-rays is dominated by the upscattered thermal disk 
emission but, because of the low densities, lies under the 
data and also has the wrong shape. In this and all other 
figures in this paper, we do not include IC of photons 
from far out in the jet, since they are orders of magnitude 
under the data and will not affect the fit. Changing the 
disk geometry so that the thin disk underlies the base of 
the jet/corona may provide more IC emission, but would 
still not address the fundamental problem of explaining 
the radio thru IR data with reasonable parameters. 

In summary, the energy budget in the non-accelerated 
case is not large enough to create enough synchrotron to 
fit the radio for a non-relativistic plasma, since the re- 
quired Qj for the already weakly relativistic temperature 
of Fig. |l|a is super-Eddington. Considering that the radio 
through IR bands are likely due to jet synchrotron in the 
LHS, we conclude, based on all the above difficulties, that 
a model which does not include particle acceleration in 
the jet is unlikely. The extent to which the energy budget 
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Fig. 3. The same model for the 1981 data set presented in 
Fig. ^3, in units of vly to show the distribution of power 
across the frequency bands. Above lO^'* Hz one can see an 
analogy to the "camel back" feature known from blazar 
spectra. 

falls short is in part due to assumptions of the equipar- 
tition of energy, as described earlier. However, corona or 
other accretion flow models which attempt to explain the 
low-frequency data without a standard jet solution do not 
seem very feasible. 

4.3. 1981 data set, with acceleration 

Therefore we want to focus on solutions where the parti- 
cles in the jet are accelerated into a power-law distribu- 
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Fig. 2. Four fits for tlie 1981 data set shown in Fig. |l|, with the radio point extrapolated from the correlation curves 
of COO and C02. The solid thick line is the total spectrum, the dotted line is the multi-temperature blackbody outer 
disk plus single blackbody irradiation contribution, the short-dashed is the synchrotron emission from the jet before 
the shock acceleration region, the dot-dashed line is the synchrotron emission after acceleration, and the long-dashed 
line is the inverse Compton upscattered disk and jet photons, (a) The fit for the mildly relativistic case of Te = 2 x 10® 
K, which requires an unrealistically large jet power of Qj ~ 0.6iEdd- (b) A fit with the highest temperature electrons 
allowed which can give a good fit to the broadband spectrum. To = 7 x 10®. This solution also gives a more realistic 
jet power of Qj O.SLEdd- (c) The case for the same temperature but with a much smaller inclination angle, requiring 
only Qj ~ 3 X 10~^LEdd, but which cannot fit the extrapolated radio point, (d) The resulting fit from a power-law of 
particles at the base of the jet, rather than a thermal distribution, as may be expected near an accretion shock. This 
solution also requires a large jet power, Qj ^ O.SLEdd: however. 



tion, as is typical in AGN. This solves the energy budget 
problem for the jet power, although we still need a very 
luminous thermal contribution. 

We show four representative fits in Fig. ^, all with a 
nozzle radius of 3rg, which is more appropriate for the 
assumption that the jet originates at the edge of an ac- 
cretion fiow. As discussed above, we fix = O.li^Edd 
with Tin — 10^ and Tout — IQ^rg and then fit the opti- 
cal data with the irradiated component, leaving the in- 
tercepted fraction of Qj, Tbb, as a free parameter in the 
fit. 

In panel (a) of Fig. |^, we show the resulting spectrum 
from a jet with mildly relativistic electrons, similar to the 
model shown in Fig Ilia. The short-dashed line shows their 



contribution before encountering the acceleration region. 
Compared to Fig. ^ the nozzle "hump" lies at a higher 
frequency due to the more compact nozzle scale, and the 
slope is more inverted primarily to the greatly increased 
angle to the line-of-sight necessitated by the radio/IR fit, 
9i = 55°. The slope is fiatter in the post-accelerated com- 
ponent (shown with the dot-dashed line), because of the 
effects of reacceleration and particle losses. Synchrotron 
from the accelerated particles dominates the entire spec- 
trum with the exception of the optical, with the IR and 
X-rays stemming from the region near the start of the ac- 
celeration Zsh, while the radio comes from much further 
out along the jet. For a fixed To, Qj and 6i are found by 
fitting the radio-IR spectrum, which then uniquely deter- 
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mines Zgh by fitting the IR turnover and extension to the 
X-rays. This also determines the slope of the accelerated 
particles, p. The jet parameter values for this mildly rela- 
tivistic case are: Te = 2 x lO^i^, Qj = 0.62LEdd, Oi = 55°, 
Zgh = 2.15 X and p — 2.15. The IC emission is dom- 

inated by upscattered thermal disk photons and is well 
under the synchrotron contribution. 

The problem with this model is that — as in the model 
from Fig. |^a — the required power in the jet Qj is unreal- 
istically large, in order to compensate for the low electron 
temperature. This imbalance also leads to the magnetic 
field energy density in the jet being a few orders of magni- 
tude super-equipartition with the hot, radiating electrons. 
While many jet formation models assume the magnetic 
energy density is in equipartition with the kinetic energy 
near the base of the jet (e.g., Blandford fc Payne 1982 ; 
Li et alj |1992| ), which for a maximal jet would lead to 
magnetic domination, this may be rather extreme. 

In Fig. ||b, we consider the highest temperature elec- 
trons allowed by this model, above which the X-ray pre- 
dictions (from either synchrotron or IC from the jet) are 
too high, Tc = 7 X 10^ K. The jet free parameters are 
now: Qj = 0.27LEdd, Oi = 55°, z^h 1-75 x lO^r^ and 
p — 2.15. Raising the temperature thus serves mainly to 
reduce the power requirements and changes the location 
of the shock region slightly. With the same SD parameters 
as in (a), the irradiation component now comprises 0.4Qj, 
still rather large but again given the uncertainties in the 
mass and outer radius, more viable. The hotter electrons 
now also contribute to a greater degree via IC upscatter- 
ing, contributing some of the soft X-ray emission. 

While there is no measurement of the inclination an- 
gle, some papers favor a lower value than what we find for 
our best fit above (e.g., Cowley et al. 2002| ). In Fig. |^c, 
we show what our model predicts for Qi — 15° and the 
same temperature as in panel (b). The decreased angle 
greatly reduces the power requirements in the jet, but at 
the same time flattens the predicted radio-IR slope signif- 
icantly, making a good radio/IR fit impossible. The pa- 
rameters are Qj = 3.2 x 10~^LEddj ^sh = 2.15 x Vd^Vg and 
p — 2.15. Although this radio point is only from an ex- 
trapolation of the radio/X-ray correlation relation (COO; 
C02), the slope is remarkably well-defined at these higher 
luminosities. On the other hand, Bi — 55° does give good 
fits to most of the other data sets discussed below. Because 
the beaming reduces the power requirements for the jet, a 
low inclination angle also means that more power must ir- 
radiate the outer disk than is available, making this model 
problematic. 

Finally in Fig. ||d we consider the case of a power-law 
of particles entering the base of the jet from the accre- 
tion flow instead of thermal particles, as may occur if ac- 
celeration also takes place at an accretion shock. Instead 
of defining a temperature, we take the minimum electron 
Lorentz factor 7o,min = 1, leaving 7c, max a free parame- 
ter. In order to prevent the nozzle emission from violating 
the constraints of the soft X-ray data, we find 7c, max ^ 10, 
which would require extreme radiative cooling to truncate 



the accelerated power-law for this scenario to be realistic. 
This type of scenario may be necessary for certain low- 
luminosity AGN (see, e.g.. Yuan et al. 2002). The other 
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jet parameters are Qj = 0.55LEdd, 2sh — 2.15 x 10"r 
di — 55° and p — 2.15. The high required jet power means 
that the irradiated component needs only 0.02(5j to ac- 
count for the optical data. 

In this last model, the IC emission still lies under the 
synchrotron component, but not by much, and has the cor- 
rect slope. In a different geometry where the SD underlies 
the corona/jet base, this may provide the best alternative 
for a solution where the IC dominates the synchrotron 
emission in the X-rays. 

In all four model runs shown here, we needed a rather 
large fraction of the jet power to be intercepted by the 
outer disk and reradiated. If we allowed the SD to radi- 
ate 0.2LEdd away, we could lower the necessary Lbb, but 
then again we would be basically at the Eddington accre- 
tion rate, as discussed for the non-acceleration case. This 
problem is clearly not very constrained for this source, 
and a central mass of ^ 15 instead of 5 Mq (as in, e.g., 



GRS 1915-H05; |Greiner et al.|p001| ), would easily resolve 
these issues. 

In Fig. m we show the same model as in Fig. ||b with 
the flux multiplied by the frequency to show the distribu- 
tion of power across the wavebands. One can see that most 
power seems to fall in the disk component, with the syn- 
chrotron then creating the second peak around the cutoff. 
The IC component is in fact dominated by the upscattered 
disk photons, but this is hidden beneath the synchrotron 
emission. The self-Comptonized synchrotron component 
peaks at 10^'^ Hz, or ^ 400 MeV, corresponding to the 
synchrotron peak upscattered by the hot electrons near 
the shock. At high frequencies, the shape is similar to the 



cam el back" spectra seen in blazars (e.g., [Fossati et al 



1998), and in fact at this highest observed flux level is 



above the sensitivity of GLAST for a one year all-sky sur- 
vey at 100 MeV. The less powerful epochs of this source 
presented in the next section may not be so observable, 
however, as they are orders of magnitude lower. 

In conclusion, we have shown that in order to ex- 
plain the low- frequency data from the 1981 observation 
of GX 339—4 , a jet solution must include some form of 
particle acceleration. And further, once that condition is 
fulfilled. X-rays from jet synchrotron (in addition to IC) 
are unavoidable as long as IC emission is not the main 
cooling channel. But under the assumption of a disk ge- 
ometry where a SD transforms to an optically thin, hot 
accretion fiow at some inner radius, synchrotron cooling 
will dominate in a maximal jet. It is possible, however, 
that in the alternative geometry of a corona/jet base lying 
directly on top of the SD could lead to so much IC cooling 
that the synchrotron cutoff would occur below ^ 10^^ Hz, 
thus accommodating a dominant IC contribution in the 
X-rays. This would not, however, explain the interesting 
coincidence of IR turnover/X-ray spectrum seen here and 
now also in several sources in addition to XTE J1118-I-480. 
It may also be difficult to explain the radio/X-ray corre- 
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lation with an IC model (see Sect. ||). We will, however, 
soon explore combined jet/corona solutions using lower 
temperature cyclo-synchrotron and IC scattering in dif- 
ferent geometries. 

We have explored the relativistic electron parameter 
range for the 1981 data set and find that the most rea- 
sonable parameters are obtained for the case of a fully 
relativistic initial thermal electron distribution. Besides 
requiring too much power, a mildly relativistic initial elec- 
tron distribution would also result in the magnetic energy 
density being much higher (and further out of equipar- 
tition with) the accelerated particles at the shock, more 
so than would likely be expected for magnetic confine- 
ment. For these reasons, in this paper we consider only 
the higher temperature solutions for all fits, correspond- 
ing to the type of fit shown in Fig. ||b. 

4.4. Other data sets 

For the 12 data sets compiled during the other LHS 
episodes of GX 339—4 , we can either optimize each fit 
individually leaving all parameters free, or fix as many pa- 
rameters as possible to see if essentially the same model 
can account for all observed spectra. We here choose to 
take this latter approach, in order to explore the funda- 
mental physics of the model. We will fix some of the pa- 
rameters to those best constrained by the 1981 set and 
then study the affect of changing only two jet parame- 
ters. Our emphasis is on understanding the radio/X-ray 
correlations, and thus the jet component. However, we ad- 
dress the SD contribution as well, and attempt to test the 
jet/disk connection by finding a scaling solution. We set 
-^BB = O.lQj for all fits, and then fit id to the optical 
data, when present. 

While it seems clear that there are small variations in 
the soft X-ray spectra of these data sets (this is discussed 
in NWD, who fit a broken power law — we are fitting the 
lower energy segment), we choose to fix the accelerated 
spectral index to p = 2.15 as in 1981. Similarly, while there 
are obvious changes in the SD contribution, they seem con- 
sistent with the likely underlying accretion changes that 
also affect the jet since both components decrease in lu- 
minosity with increasing time. We keep Tjn fixed at 10^ 
K and Te fixed at 7 x 10^ K based on the 1981 set, as a 
starting point. We will assume that the orientation of the 
jet remains fixed at 55° (statistically 57° is the most likely 
angle) and that the base of the jet is fixed at ro = 3rg. 
Thus the only jet parameters allowed to vary are Qj and 
Zsh, and the disk parameter is then adjusted to fit the 
thermal data, when present. 

Figs, ^a & show the broadband spectral fits to 
the data, with the corresponding fit parameters given in 
Table ^. For the data sets without optical data constrain- 
ing the disk, is rather arbitrary, and we scaled the 
luminosity to the jet power. 

Although the fits are surprisingly good given the lim- 
ited fitting parameters, there are obvious features in the 




Fig. 5. Correlation between the input jet power Qj and 
the location of the first acceleration zone Zsh for the fits in 
Figs. H and ^. The dashed line is to guide the eye for the 
upward trend in the data sets. The 1981 point represents 
the fit shown in Fig. |^b. The 998 data set is not strictly 
simultaneous, and the radio data are simultaneous with 
a small reflare in the hard X-rays, which may explain its 
excess. 99/ may also be affected by this flare. 

X-ray waveband such as a flattening at higher frequencies 
and peaks which our model cannot account for. These fea- 
tures are likely due to the impacting X-rays on the cooler 
material of the SD, but we have not yet incorporated this 
into our model. The softest X-rays should not be affected 
by these additional components, however, so we focus our 
fitting at the lower X-ray frequencies. For the lowest fiux 
level data sets of 1999, which we label 994 — 998, as well as 
2000, the error bars are large enough to make the fitting 
more difficult. 

The most important feature to note is that the 
radio/X-ray correlations detected in this source (COO; 
C02) can be explained by the variation of only one main 
parameter in the jet model: the input power. If we plot 

against the jet power Qj (Fig. ||), we see that there is 
a slight dependence (zsh (5"°'^'^^), but this likely results 
from our rather artificially freezing of the other parame- 
ters. This small dependence is rather amazing when one 
thinks that Qj is changing by about three orders of mag- 
nitude over the same range, if 1981 is included. 

5. Understanding the radio/X-ray correlations 

It seems likely that the real source of the variations is the 
power fed into the jet. In Fig. ^ we plot the prediction 
of our model against the data for the correlation between 
the 8.6 GHz radio fiux and the integrated 3-9 keV X-rays 
(C02). We fix ~ 1-8 x lO^r^, to match the highest 
flux data point of 1981, which is not shown on this flg- 
ure because the radio was not simultaneously measured. 
Then we adjust only Qj while holding every other param- 
eter constant, in order to fit the radio flux, and plot the 
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Table 2. Fit Parameters for Figures 3a and 3b 
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resulting integrated 3-9 keV X-ray flux. In removing the 
Zsh oc Qr°'^'^^ dependence mentioned above, we are likely 
compromising the quality of the fit to some degree, but 
this could be compensated for if we allowed particularly 
the spectral index to vary. 

By just changing the jet power we obtain a very good 
fit to the correlation data, with the exception of the two 
data sets 997 and 998. In fact, the radio data in 998 was 
taken a few days after the X-rays, and coincided with a 
small reflare in the hard X-rays (see Fig. 15 in COO). This 
is the likely explanation of the excess radio flux. Similarly, 
the reflare may also have an affect on 997, if there is, e.g., 
a lag between the changing fluxes. 

The results shown in Fig. ^ provide strong support for 
a jet synchrotron model, since it correctly predicts the 
radio/X-ray correlation over orders of magnitude in flux 
levels, almost exactly matching the data. The fit is better 
at the higher flux end, where the error bars are smaller. 
The exact dependence (excluding 997 and 998) goes as 
Fx oc F^, where m = 1.41, Fr is the 8.6 GHz radio flux 
and Fx is the integrated 3-9 keV X-ray flux. 

The value for m in fact follows directly from ana- 
lytic predictions of the jet model. 



Falcke & 



Eq. (52) ^ 

Bierm£,nn (1995) shows|^ that if all parameters except the 
power are fixed, as we have done to produce Fig. ^, the 
frequency where the optically thick fiat-to-inverted spec- 
trum turns over to the optically thin regime depends on 
the power as 

i^SSA oc Q^^^. (3) 

Similarly Eq. (56) (ibid.) shows that under the same freez- 
ing of other parameters, the flux at the turnover depends 
on the power as 



FssA oc Q 



17/12 



(4) 



These relations can then be normalized using the 1981 
data set, where the flux and turnover frequency are visible. 



In Falcke i 



Biermann (1995), the quantities which we label 

qj/iLie and L^^^obs, 



fssA, Qi and Fssa correspond to !/s,obs 
respectively. 



Knowing these relations, and the spectral indices of 
the optically thick and thin components from the data, 
one can roughly calculate the expected correlation slope. 
Because the spectral indices for the radio-IR, cirir, and 
X-ray, ax, do vary slightly over GX 339—4 's observation 
history, taking constant values will make this less exact. 
We now want to calculate the slope between two points 
along the correlation curve, analogous to the model plot- 
ted in Fig. ^. Each point represents a different value of 
Qjct, corresponding to different values of v^sA and Fssa- 
The 8.6 GHz radio flux is then 



Fs.e — 



^^^^ ("8.6 GHzV 
and then using Eqs. ^ & 



(5) 



logio ^8.6 



17 2 
logio + Y2 ^°^io ^ 3"RiR logio Qj, (6) 



where Ci absorbs the exact dependences of the flux and 
frequency at the turnover on Qj, and other constants. 

Similarly, the X-ray flux at any frequency (as long as 
it falls on the power-law) is 

X QX 



Fx — F' 



SSA 



(7) 



which after using the forms from Eqs. (|^) & (Q) and inte- 
grating from 3-9 keV becomes 

f 17 

logio / Fxdiy = logio '^2 + — log^ Qj- 

2 

3 ax logio Qj, (8) 

where again C2 includes the constants of integration and 
the dependencies on Qj. 

For any two values of the jet power, Qa and Qb, the 
correlation slope is then 

logio / -Px,B - logio / -Px.A 



logio Fs.G.B - logio ^8.6,A 
17 2, 



12 3 



ax 



17 2„, ' 

12 - aaRiR 



(9) 

(10) 
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log,„ 8.6 GHz Radio Flux (mjy) 

Fig. 6. The model-predicted radio (8.6 GHz)/X-ray (in- 
tegrated 3-9 keV emission) correlation (solid line). The 
data are from COO & C02. For this figure, Zsh is fixed at 
~ 1.8 X lO^Tg to match the 1981 highest flux result, and 
then the only changing parameter is the power input into 
the jet Qy The arrow represents 3 — cr upper limits. 

after algebraic cancellations. Taking the values from the 
power-law fits to the 1981 data set from CF02, aRiR ^ 
0.15 and ax ^ —0.6, we find acor = 1.38, which consider- 
ing the simplifications in the model and our ignoring of the 
spectral variations, is surprisingly close to the value of 1.4 
found for GX 339—4 . In addition, this type of radio/X-ray 
correlation has also been seen in V404 Cygni, with exactly 
the same slope ( pallo et al1|2002| ; Gallo et al., MNRAS, in 
prep.). Our preliminary models of this source also show a 
turnover in the IR band and have similar jet parameters 
(Markoff et al. 2002, in prep.). 

The fact that the well-determined slope m extrapolates 
to the lower flux data also gives us greater confldence in 
the subtraction of the putative "background source" dis- 
cussed in Sect. 2. Note also that the lowest flux point is 
from BeppoSAX, and therefore does not rely on the sub- 
traction of a background source. 

It is possible that an IC corona model could also ex- 
plain the correlation slope, if the radio flux scales with 
the power (assumed to be linearly related to M) as ex- 
pected from the jet and the X-ray flux scales with AP. 
This may be the case in some optically thin accretion solu- 
tions. However, it has yet to be shown that a self consistent 
solution can be found which also well flts all the broad- 
band data available for this source. Similarly this would 
not explain why the X-rays trace back to the turnover in 
the IR as 1981 data, or why most if not all LHS sources 
show this turnover coincidence. 



6. Conclusions 

We show that a model comprising a dominantly syn- 
chrotron jet component, in combination with an optically 
thin accretion flow transitioning to a standard thin disk. 



is able to explain the broadband spectral data from 13 
observations of GX 339—4 . At the same time, the model 
can easily explain the m ^ 1.4 slope of the radio/X-ray 
correlations by changing only the power input into the jet. 
This input power, Qj , is assumed to be proportional to the 
accretion rate, and varies from LEdd 0.003 — 0.3 for a 
5Mq black hole. Above the highest luminosity, the source 
would likely transition to the HSS. 

This model assumes an underlying disk geometry in 
which the standard thin disk exists beyond ~ 100 — 
1000 Tg, and has a relatively low accretion rate. In a sit- 
uation where the SD extends all the way down to the 
last stable orbit, the photon fleld at the base of the jet 
may be high enough to affect our conclusions, which are 
dependent on synchrotron losses dominating the cooling. 
If the available photon fleld for inverse Compton upscat- 
tering becomes large enough, it may truncate the accel- 
erated particle distribution and lessen the extent of the 
synchrotron emission. 

With the exception of this model, the non-thermal X- 
ray spectral component in all XRB states has been mod- 
eled in terms of inverse Compton processes in a hot corona, 
located either above or within the standard thin disk. 
These classes of models have significant success explain- 
ing both the spectral features and timing characteristics 
in the X-ray waveband, but fail to address the radio/IR 
spectra which in at least GX 339—4 , V404 Cyg and Cyg 
X-1 seem intimately linked to the higher frequencies. 

However, the existence of a hot, magnetized plasma at 
the base of the jet leads to obvious considerations of a uni- 
fication scheme between it and this concept of a corona, 
which until now has not been explicitly observed. Both 
scenarios (jet and corona) address important features from 
the observations, but either cannot or have not yet at- 
tempted to account for everything. We feel that a con- 
catenation of these two components would be fruitful as 
the next stage of our investigation. It is, however, impor- 
tant to note that under reasonable physical assumptions, 
the jet can easily produce X-ray emission via synchrotron 
radiation, and has significant success explaining the data 
when it does. This possibility must therefore be considered 
in spectral modeling; whether it really dominates the en- 
tire spectrum is a question of the environment and local 
acceleration conditions. It may be that it indeed domi- 
nates only in those sources which show only weak reflec- 
tion features (e.g., XTE J1118+480; [Miller et al.]^00l| ), 
but it likely has a wider impact in all sources than has so 
far been appreciated. 

Interestingly, the parameters we derive for the power 
and location of the shock are roughly consistent with ear- 
lier results modeling XTE J1118-f480, which because of 
the higher temperature (Te = 2 x 10"'^° K) resulted in a 
lower necessary power input into the jet. In this source we 



found 



10^ 



as compared to the ~ 10 found 



here, however considering that the jet can extend be- 
yond lO^^Tg, this range is quite small. In our recent mod- 
eling of other LHS BHC sources, we are in fact flnd- 
ing that all sources with simultaneous radio/X-ray data 
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seem to require acceleration to begin in this same range, 
and we will discuss the physical implications of this else- 
where. For these sources wc also find, as in GX 339—4 and 
XTE J1118+480, that if the X-rays are traced back to 
lower frequencies, the optically thick-to-thin turnover al- 
ways occurs in the IR range. If this turnover coincidence 
is real, it offers us the chance to explore the physical con- 
ditions at the acceleration zone of XRB jets. 

This work illustrates how critical simultaneous multi- 
wavelength observations are to increasing our understand- 
ing of the physics in these sources. With the advent of 
higher energy missions such as GLAST, INTEGRAL and 
ASTRO E-2, soon we can hopefully begin to probe the 
contribution of these smaller jets to the hard X-ray and 
7-ray bands. 
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Fig. 4. (a) Composite figure with 6 data sets with labels corresponding to the data sets of Table ^ taken from W99, 
COO, C02 and NWD. The solid thicker line is the total spectrum, the dotted line is the multi-temperature blackbody 
representing the outer thermal disk plus a single blackbody for irradiation, the short dashed line is the synchrotron 
emission from the jet before encountering the shock acceleration region, and the long-dashed line is that from after, 
and the dot-dashed line is the IC jet and disk emission upscattered by the jet plasma. Fit parameters are shown in 
Table H 



16 



Markoff et al: The Role of Jets in GX 339-4 




10 15 20 10 15 20 

logio V (Hz) logio V (Hz) 



Fig. 4. (b) Composite figure with 6 data sets with labels corresponding to the data sets of Table ^ taken from W99, 
COO, C02 and NWD. The arrows represent 3 — cr detection upper limits. The lines are the same as in Fig. For the 
'off state in 2000, the optical data are from three months prior to the X-ray and radio observations, and should be 
seen as upper limits. 



